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ly-j ; Abstract 

We have studied measurements of the cross section for e + e~ — > hadrons for center- 
■ of-mass energies in the range 20-209 GeV. We find an apparent excess over the pre- 

dictions of the Standard Model across the whole range amounting to more than 4a. 
As an example, we compare the data to predictions for a light scalar down-type 



■ quark which fit the excess well. 

o 

2t ' 1 Introduction 

JJjj ■ Measurements of the inclusive cross section for e + e~ — ► hadrons have been a staple of e + e" 



Q-i 



collider experiments for many years. They have been essential for testing the predictions of 



QCD, and for measurements of the properties of the Z boson which underpin the electroweak 
sector of the Standard Model (SM). It is often said that the SM survives all confrontations 
with data, and indeed, this agreement is a very important constraint on models for new 
^ ■ physics beyond the SM. For example, one of the strengths of the Minimal Supersymmetric 
extension of the Standard Model (MSSM) is the fact that it does not disturb this agreement 
though its parameters may be varied over wide ranges. 

It has been noted several times that the Lep 2 measurements of e + e~ — > hadrons tend 
to exceed the predictions of the SM P . It is also well known that the hadronic cross section 
at the Z peak, cx° ad , slightly exceeds the SM fit |21E]- Neither of these excesses is statistically 
significant. 

We have examined the data on e + e~ — > hadrons from experiments that ran below the 
Z peak using a compilation amassed by Zenin et ai, [4J, and find that these measurements, 
taken together, also are in excess of the SM prediction. We developed a likelihood method to 
quantify the significance of the excess in four ranges of the e + e~ center-of-mass energy, y/s, 
and then combined the likelihoods. The net excess has a significance of over 6a, if theoretical 
uncertainties and correlations for experimental systematic uncertainties are ignored. When 
a conservative theoretical uncertainty is applied, the significance reduces to 4.4cx, and when 
the correlations among experimental measurements are taken into account, the significance 
is 4.3o", corresponding to a probability of 10~ 5 . 



2 Apparent Excess in e + e — > hadrons 

We compare the data to three straw models of new physics: 1) a purely phenomenological 
ansatz with <r NP oc 1/s, 2) light bottom squarks (b\) with little or no coupling to the Z, and 
3) an additional neutral boson resonance (Z') with a mass well beyond the energies of the 
data. We find that the data are consistent with a light b±, and do not favor an additional Z'. 

The structure of this report is as follows. First, we discuss the data used in this study, 
and then the SM prediction for e + e~ — ► hadrons. Next we compare the measurements to 
the SM prediction, which is based on the likelihood method. After establishing the excess on 
the basis of a very simple expression for any possible excess, we compare to the expectations 
for a light bi, and also to the expectation for a heavy Z'. Finally, we draw our conclusions. 

2 The Data 

We base this analysis on data from the LepEWWG [2] and the compilation of measurements 
by Zenin et al. [I]. 

The results from the analysis of Lep data have been stable for some years. They fall 
naturally into two groups: measurements around the Z peak (88 < < 93 GeV) and well 
above (130 < y/s < 210 GeV). We rely on the reports in Refs. [21 Bj, both for the reduction 
of the measurements and for the SM predictions. A summary of the data is given in Tabled 

The data below the Z peak come from a number of experiments running at several e + e~ 
colliders. The main interest in these data comes from the need for improved estimates of 
«qed(^1) and corrections to (g — 2) M . We consider two subsets, namely, 20 < yjs < 40 GeV, 
for which Z-exchange should be unimportant, and 40 < ^/s < 70 GeV, which will show the 
onset of Z-exchange. The combined measurements in these two regions carry approximately 
equal precision. A summary of the relevant experiments is given in Table El 

These considerations lead us to designate four regions based on y/s: 

1 20 < v/i < 40 GeV 

2 40 < v/i < 70 GeV 

3 88 < v/i < 93 GeV 

4 130 < v/i < 210 GeV 

We do not consider the data at y/s < 20 GeV as some of these measurements are not very 
precise, the theoretical uncertainty from as becomes large, and the prediction for any new 
particle might be complicated due to threshold effects. Plots of the data from regions 1, 
2 & 4 are displayed in Figs. [H El and El 

There are 130 measurements for regions 1 & 2. The statistical precision of the typical 
measurement is ~ 5%, and the total systematic uncertainties are roughly the same as the 
statistical errors. In some cases the systematic uncertainties are not well described. Gener- 
ally speaking, they involve the luminosity measurements, the hadronization model, and the 
trigger efficiency 
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Figure 1: Data for region 1 (20 GeV < \/s < 40 GeV) and the Standard Model prediction 
from Zfitter. 
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Figure 2: Data for region 2 (40 GeV < y/s < 75 GeV) and the Standard Model prediction 
from Zfitter. 
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Figure 3: Ratio of the Lep 2 measurements to the SM prediction from the LepEWWG 0/. 



energy 


combined measurement 


SM prediction 


difference 


deviation 


91.187 


41540 ± 37 


41478 


62 


1.68 


130 


82.1 ±2.2 


82.8 


-0.7 


-0.32 


136 


66.7 ±2.0 


66.6 


0.1 


0.05 


161 


37.0 ± 1.1 


35.2 


1.8 


1.64 


172 


29.23 ±0.99 


28.74 


0.49 


0.49 


183 


24.59 ± 0.42 


24.20 


0.39 


0.92 


189 


22.47 ±0.24 


22.16 


0.31 


1.29 


192 


22.05 ±0.53 


21.24 


0.81 


1.53 


196 


20.53 ±0.34 


20.13 


0.40 


1.18 


200 


19.25 ±0.32 


19.09 


0.16 


0.50 


202 


19.07 ±0.44 


18.57 


0.50 


1.13 


205 


18.17 ±0.31 


17.81 


0.36 


1.16 


207 


17.49 ±0.26 


17.42 


0.07 


0.27 



Table 1: Summary of the Lep measurements, and SM predictions, from Ref. 0/. The com- 
bined measurements and SM predictions are given in pb. The 'difference' is (measurement - 
prediction), and the 'deviation' is the difference divided by the experimental error. 



experiment 



(GeV) 



publication 



CELLO / PETRA 



34.9, 42.7 

22.0, 33.8, 38.3, 41.5, 43.5, 44.2, 46.0, 46.6 



Phys.Lett 144B (1984) 297 
Phys.Lett 183B (1987) 400 



JADE / PETRA 



22.0, 25.0, 27.7, 29.9, 30.4, 31.3, 33.9, 34.5, 

35.0, 35.4, 36.4, 40.3, 41.2, 42.6, 43.5, 44.4, 45.6, 46.5 
22.0, 27.6, 30.8 



Phys.Rep. 148 (1987) 67 
Phys.Rep. 83 (1981) 151 



Mark-J / PETRA 



31.6 

22.0, 27.6, 30.0, 30.7, 31.6, 33.0, 34.0, 35.3, 35.8 
34.8 

22.0, 25.0, 
40.3, 41 



30.6, 33.8, 34.6, 35.1, 36.3, 37.4, 38.4, 
5, 42.5, 43.5, 44.2, 45.5, 47.5 



Phys.Lett. 85B (1979) 463 
Phys.Rep. 63 (1980) 337 
Phys.Lett. 108B (1982) 63 

Phys.Rev. D34 (1986) 681 



TASSO / PETRA 



22.0, 27.7, 
22.0, 25.0, 
22.0, 25.0, 
41.4, 44.2 
22.0, 35.0, 



30.9, 31.6 

27.6, 30.2, 31.0, 33.0, 34.0, 35.0, 36.0 

27.7, 30.1, 31.5, 33.5, 34.5, 35.5, 36.7, 43.1 

43.7 



Z.Phys. C4 (1979) 87 
Phys.Lett. 113B (1982) 499 
Z.Phys. C22 (1984) 307 
Phys.Lett. 138B (1984) 441 
Z.Phys. C47 (1990) 187 



Mark-II / PEP 



29.0 



Phys.Rev. D43 (1990) 34 



MAC / PEP 



29.0 



Phys.Rev. D31 (1985) 1537 



AMY / TRISTAN 



50.0, 52.0, 54.0, 55.0, 56.0, 56.5, 57.0, 58.5, 
59.0, 59.1, 60.0, 60.8, 61.4 



Phys.Rev. D42 (1990) 1339 



TOPAZ / TRISTAN 



50.0, 52.0, 
59.1, 60 
57.4, 58.0, 
57.8 



54.0, 55.0, 56.0, 56.5, 57.0, 58.3, 
0, 60.8, 61.4 

58.2, 58.5, 58.7, 59.0, 59.2, 59.5, 59. 



Phys.Lett. 
Phys.Lett. 
Phvs.Lett. 



234B (1990) 525 
B304 (1993) 373 
B347 (1995) 171 



VENUS / TRISTAN 



50.0, 52.0 
63.6, 64.0 



Phys.Lett. 
Phys.Lett. 



198B (1987) 570 
B246 (1990) 297 
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3 Standard Model Predictions 



The SM predictions for the inclusive process e + e — > hadrons are relatively simple. Let M 



f 



be the mass of a particular quark, so (3 = ^/l — AMj/s is its velocity. Denote its electric 



charge by Qf, and weak isospin by // = ±|. The fundamental constants are the Fermi 
constant, Gp, the QED coupling oqed, and the strong coupling as- The latter two run as 
a function of s. Finally, the weak mixing angle relates weak and electric couplings, and can 
be defined according to 

2 _ 7ra QED 



Gf Mi 



V2 sin 2 9w cos 2 9 



w 



where Mz is the mass of the Z boson. Since aQED(s) runs, so does sin 2 9w, according to this 
definition. The vector and axial vector couplings of the quark / to the Z can be written 

v _ l(-2Q f sm 2 9 w a I_[ 
^ 2 sin 2 9w cos 2 9w * 2 sin 2 9w cos 2 9w 

and the Z propagator in the lowest order is simply [3] 

s 

s-Ml + iMzTz 



Xo(s) = - — ^ 2 nnrr-F ( 3 ) 



where V z is the Z width. 

The amplitude for e + e~ — > qq is the sum of the amplitudes for 7* and Z exchange, so 
the cross section is the sum of three terms. After integrating over all angles, 



cr(e + e — > qq) 



2 * a QED 

3s 

Q) - 2v e v f Q f Re Xo (s) + (v 2 e + a 2 ) L) + a) f^f^ ) ) IXo(s)P 




(4) 



where C$ is a QCD correction factor depending on as and the number of colors, iV^ = 3. 

Equation (0J pertains to the simple Born approximation when «qed, and hence sin 2 6 W , 
are fixed, and Cs = Nc- The improved Born approximation allows for running q;qed and 
sin 2 9w, for finite corrections to Cs/Nc, and for an energy dependent width: Tz = T° z s/Mg, 
with T° z independent of s. We list the values for all required constants in Table El 

We use version 6.36 of the Zfitter program jHj to compute the SM predictions in regions 
1 & 2. This program has been developed by several authors over many years, and is one 
of the two main programs employed by the LepEWWG. For an in-depth discussion of the 
theoretical accuracy of the predictions from Zfitter, see the paper by Bardin, Griinewald, 
and Passarino [7j. 

The routine for computing the s-dependent hadronic corrections to «qed was recently 
improved by Jegerlehner jH], who provided us with nearly-final version. We use this new 
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parameter 


value 


1/«qed(0) 
Aqcd 


137.036 
2131^ MeV 


M z 


91.1876 ±0.0021 GeV 


1 z 


2.4952 ±0.0023 GeV 


M b 


4.5 GeV 


M c 


1.6 GeV 



Table 3: Constants required for the evaluation of the cross section in the improved Born 
approximation. The given value for Aqcd leads to as(M z ) = 0.1184 ± 0.0072. The uncer- 
tainties on ttQED; Mb and M c are irrelevant for this analysis. 

routine for our SM computations. The change in Aa^ had ^ is generally less than 2%, and so 
the impact on the cross section is negligible. 

We wrote our own code to compute the cross section for e + e~ — > hadrons, based on 
Eq. For cyqed(s) we used a routine from H.Burkhardt 0, and for the QCD correction 
factor, Cs, we used an expression published by S.Bethke jTUj: 




with Ki = 1 and K2 = 1.4. The exact values of K\ and K2 depend slightly on the quark 
masses, and we consider the uncertainties of these coefficients below. 

We compared the results of our code to those of Zfitter. There is agreement at the 
level of 0.2 % after accounting for the fact that the running of the heavy quark masses is 
implemented in Zfitter but not in our code; the predictions from Zfitter are slightly 
higher than those from our code at the level of 0.1-0.2 %. Also, box diagrams contribute 
about 0.5 % to e + e~ — > hadrons above the W + W~ threshold, which is not included in the 
simple expression Eq. (j3J). For the comparison of the SM prediction to the data, we use the 
values from Zfitter, as this program is extremely well tested and reliable. Our code is 
useful for investigating sources of theoretical errors, as discussed below. 

We must assess the uncertainty on the SM prediction in order to make the comparison 
to the data. The studies by Bardin et al. [7j indicate that the theoretical accuracy is 0.2% or 
better, in regions 1 & 2. For region 4, the Two-Fermion Working Group reports a theoretical 
uncertainty of 0.26% [TT] . They also report values for 1/«qed which imply an uncertainty 
on a(e + e~ — > hadrons) of 0.14%, and we obtain similar results. 

The uncertainty coming from the inputs must also be estimated. For ttQED(s), we must 
consider both the uncertainty on its measured value at yfs ~ 0, and the uncertainty from the 
running, which is dominated by the uncertainty on Aa^ hlld \ The uncertainty from «qed(0) 
is negligible. We estimate the uncertainty from the running by modifying the coefficients in 
the Burkhardt routine to correspond to Aa^ 1 ^ = 0.02761 ±0.00036 at the Z peak [U]. This 
is a conservative estimate and is the same used by the LepEWWG in their analyses of Lep 
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data. Other estimates suggest a significantly smaller uncertainty ^2|- Furthermore, we take 
the uncertainty from «qed on the SM cross section to be fully 5(a) /a ~ 2 5(«qed)/«qed, 
ignoring the fact that uncertainties will be smaller for the terms involving Gf, as we have 
explicitly verified using our code based on Eq. (J1J. 

The uncertainty coming from the QCD correction also has two pieces: the uncertainty 
on ctg(s) and the uncertainty on the coefficient Ki- We do not consider an uncertainty on 
K\ as this would be redundant with the uncertainty on as- As discussed by S.Bethke |TU] . 
there is some spread in the values for as coming from various sets of measurements. There 
is a tendency for the lower energy measurements, such as deep-inelastic scattering, to give 
a value of as(M^) lower than that coming from the data at higher energy, mainly from 
Lep. A typical average is as(M%) = 0.1184 ± 0.0036 [TP] . The central value is more or less 
mid-way between the low values favored by DIS measurements and the high values favored 
by Lep. In the context of the light sbottom model discussed below, as will run more slowly 
resulting in a slightly higher value at s = Mf [T^]. We choose to keep the central value 
(0.1184), but double the uncertainty to ±0.0072. In the parametrization of the running of 
as by Bethke jTUj, this corresponds to the low-energy parameter Aq CD = 213^70 MeV. We 
evaluated the impact of this uncertainty on the SM cross section using our code and find it 
to be on the order of 8(a) /a ~ 0.3%. We also varied the coefficient k 2 by 10% and found 
the impact to be negligible. 

We evaluated the uncertainties coming from the other input parameters, namely, Mz, 
T° z , Mf, and M c , and found them to be completely negligible. 

A summary of the theoretical uncertainties for the entire range of yfs is given in Fig. |3J 
The dominant error comes from the running of as, which has been very conservatively 
estimated. We display also the linear sum of the uncertainties, and use this as the guide 
for setting the theoretical uncertainty. For the range 20 GeV < ^fs < 75 GeV, we take an 
uncertainty on the order of 0.5%, which falls slightly from 20 GeV to 80 GeV, as indicated by 
the heavy line. For region 3, we do not include a theoretical uncertainty, as this parameter 
is a fit result, and the reported theoretical uncertainty is exceedingly small, amounting to 
0.012%, or 5 pb [7j. In region 4, we take an uncertainty of 0.55%, as indicated by the second 
heavy line. 

We emphasize that these theoretical uncertainties are conservative, and one could argue 
that more agressive estimates closer to 0.3 % would be justified. 



4 Comparison of the Data to the SM Prediction 

We develop a series of comparisons, starting with a simple weighted average of differences and 
ending with a likelihood analysis which incorporates theoretical and experimental systematic 
uncertainties, and correlations among errors. 
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4.1 Mean Deviations 

We compare the SM prediction to the data in each of the four kinematic regions defined above. 
First, we compute a simple x 2 an d mean deviation of the data from the SM prediction. Let 
the measurements be yi with uncertainties rji, and the SM prediction, y SM , which varies 
with s. Then 



X 



and the mean deviation and rms are 




y 



SM- 



(6) 



A 



f ( vi - v m 




and 




1-1/2 



(7) 



At this point we ignore correlations among the measurements. We summarize the values 
for x 2 , A and a& in Table EJ No value is given for region 3 as there is only one data point 
there. Note that the mean deviation is significantly non-zero and positive, even though 
the x 2 is fi ne - This is a key point and a somewhat atypical situation. There are many 
measurements, each of which is statistically consistent with the SM prediction, and so the 
value of x 2 is good. But the x 2 t es t does not distinguish between measurements which are 
higher than the prediction and those which are lower. In the present case, it turns out that 
there are many more which are higher than lower. The mean deviation, A, is designed to 
make this evident. To demonstrate that these offsets are plausible, we form local averages 
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region 


N data points 


x 2 


A 




A/a A 


1 


67 


50.4 


6.7 


2.5 


2.9 


2 


63 


55.1 


7.9 


1.4 


5.9 


4 


12 


12.2 


0.32 


0.11 


2.8 



Table 4: Values for x 2 , A and <7 A / or regions 1, 2 & 4- A and a& are quoted in pb. 



of the data for regions 1 & 2, and plot the ratio of the rebinned data to the SM prediction 
- see Fig. The likelihood analysis discussed in the next section provides a more physical 
basis for analyzing the apparent excess indicated by the values in Table HI 



4.2 Simple Likelihood Analysis 

We wish to quantify better any disagreement of the SM prediction with the data. To this 
end, we define a likelihood function, based on the measurements yi with uncertainties rji and 
the theoretical prediction, y SM + ay NP (SM + new physics), where a is a free parameter, 

v = ft -i^- ex P {- A(w - (y SM + *y NP )) 2 ) • (8) 



i=i 



2,TT Tji V 2 Vi 



It is understood that both y and ?/ NP are functions of s. It is more convenient to work 
with the negative log-likelihood function, T = — In V . We will examine the change in T as 
a function of a. 

The quantity a y NP represents a contribution from new physics. As the simplest ansatz, 
we take 

,np (s) s (10 pb) x fiOGe^ (9) 

s 

which has the generic 1/s dependence of the pair production of light particles, and we 
normalize to 10 pb at y/s = 30 GeV for convenience later. In this case the influence of Z- 
exchange is ignored. Minimizing the negative log-likelihood function J'asa function of a, we 
obtain the results (abest) summarized in Table 03 One sees that a positive contribution a > 
is favored in all cases - this comes about because the mean deviations, A are all positive 
(see Table EJ). One also sees that ignoring Z-exchange does not appear to be justified, since 
«best is much larger for region 3 than for the others. 

We compare the best values of T to those obtained for the SM alone (i.e., with a = 0) 
and see an improvement when a nonzero "new physics" contribution is included. This 
improvement can be described by a number of standard deviations (S.D.) according to 



SD 



^■(0) - ^(aw) (10) 
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Figure 5: Ratio of the measurements to the SM prediction, rebinned, for regions 1 & 2. The 
shaded band centered on 1 represents the theoretical uncertainty. 
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region 


«bcst 


-^(ttbest) 


^(0) 


S.D. 


1 


0.74 


21.89 


25.23 


2.58 


2 


1.80 


25.52 


37.50 


4.90 


3 


57.3 





1.404 


1.68 


4 


1.37 


2.00 


6.12 


2.87 


net 


1.183 


53.446 


70.255 


5.80 



Table 5: Results of the minimization of the negative log-likelihood function, T , as a function 
of a, for the simple ansatz in Eq. (0J). 



as discussed in the RPP, for example . We checked that the minimization of J- and the 
interpretation of the error on a corresponded exactly to what is obtained with a standard 
X 2 minimization, region by region. 

The values listed in Table do suggest that there is an apparent excess of data over 
the SM prediction, in all four regions of y/s. It makes sense to form a combined likelihood 
by taking the product of the likelihoods for the four regions. The corresponding negative 
log-likelihood function is then minimized with respect to the scale a, with the result: 

a net = 1.18 ± 0.20 or, y NP (s) = (11.8 ± 2.0) pb x (3 ° GeV)2 . (11) 

s 

The difference in T is 16.81, which corresponds to 5.8 standard deviations, according to 
Eq. PJ. 

The value Eq. (|llj) is significantly different from zero, but still is very small. It corre- 
sponds to only 3% of the hadronic cross section at y/s = 30 GeV. At this energy, the cross 
section for the production of a lepton pair is about 110 pb, and for a charge —1/3 quark, 
about 40 pb. Interestingly, the cross section for a scalar particle will be significantly smaller, 
due to the reduced number of spin degrees of freedom. 



4.3 Theoretical Uncertainties and Experimental Systematics 

The theoretical uncertainties depicted in Fig. 0] are not negligible compared to the deviations 
in Table EJ Furthermore, we have taken the measurements to be uncorrelated, which is not 
the case. We now rectify these deficiencies and see to what extent the significance of the 
apparent excess diminishes. 

To implement the theoretical uncertainties, we modify the likelihood function to include 
a multiplier, p, for y SM which is constrained by a Gaussian centered on one with a width 
given by the assumed theoretical uncertainty, r\ p : 



V 



N 



H ex p 



( Vi - (py SM + ay^)) 



NP\\2 



11=1 



x exp 



1 



(p 



(12) 
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where we have omitted unnecessary normalization factors. The likelihood and corresponding 
negative log-likelihood are to be viewed as functions of both the new physics scale factor, 
a, and the scale factor for the SM prediction, p. The second factor will allow a numerical 
improvement in the likelihood when a = 0, but will have no impact for the optimal abest- 
The result is a smaller difference AjF and hence a lower significance. 

We apply this new formulation of JF and find a significant reduction in AJ 7 , as expected. 
The significance is reduced from 6a to 4a. 

Next we take into account the correlation among the measurements. For the Lep 2 
data, this is straight forward, as the LepEWWG have published the covariance matrix 2J. 
For the lower energy data, no rigorous prescription is possible, so we adopt the following 
method. The data appear in twenty-one distinct publications, as listed in Table EJ It is 
reasonable to assume that the measurements reported in a given publication are correlated, 
so we take half the total systematic uncertainty to be correlated. Correlations among sets 
of measurements, however, are expected to be very small, since nearly all of the systematic 
uncertainties come from detector specific issues (such as the alignment of the luminosity 
detectors, or the measurement of the trigger efficiency), or are evaluated in different ways 
(such as the dependence on the hadronization model). 

We replace the sum over individual terms in the negative log-likelihood function by one- 
half the x 2 formed in the canonical way from the inverse of the covariance matrix. This 
X 2 function depends on a in the same way T does. The data in regions 1 & 2 are treated 
together as part of one x 2 function, and the data in region 4 is treated in a separate x 2 
function. 

Since the weight of each measurement is different once correlations are taken into account, 
the net T function must be minimized anew. Setting aside the theoretical uncertainties, the 
significance of the excess is reduced from 6a to 5a. Taking both the theoretical uncertainties 
and the correlations into account, we find a significance of 3.9a. 

5 Examples: Light Sbottoms, Extra Z Boson 

For the sake of discussion, we take the apparent excess to be genuine and compare the data 
to two ansatze for physics beyond the SM: 

1. light bottom squarks 

2. heavy Z' bosons 

The first model serves as the example of the production of new light particle in e + e~ collisions, 
while the second involves a new production mechanism of purely SM particles. They can be 
distinguished, in the absence of an analysis of the final state, by their dependence on s: the 
one falls as 1/s while the other rises approximately as s, far below the pole in the propagator. 
We do not mean to advocate either of these models, but rather use them as reasonable models 
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to test with these measurements of cr(e + e~ — > hadrons). They are both well motivated by 
other considerations. 

5.1 Light Sbottoms 

The light sbottom model comes from the proposal by E.Berger et al. to explain the excess 
of 6-hadron production at the Tevatron by postulating the existence of a light sbottom quark 
(&i), with a mass in the range 2 GeV < Mg < 6 GeV, and a light gluino (g) with a mass 
of roughly Mg ~ 12 GeV. The gluino decays only to g — ■> 6&i, and the sbottom decays via 
an .R-parity violating coupling to a pair of quarks: b\ — > ttJ, for example. For our purposes, 
the gluino plays no role - we assume that these light these bottom squarks produce hadronic 
jets through an .R-parity violating coupling (A") which would allow them to be selected in 
the measurements of a(e + e~ — > hadrons) with approximately the same efficiency as a SM 
event e + e~ — > qq. While the SM b quarks have a substantial coupling to the Z boson, the 
composition of the lightest sbottom can be tuned so make this coupling arbitrarily small. 
In the case that it is exactly zero, the cross section would have the generic s-dependence 
of Eq. (J2j). The composition is controlled by a mixing factor, cos#b , and zero Z — b\ — b\ 
coupling corresponds to cos 6*^ = 0.39. In our comparison of the light sbottom model to the 
data, we will have to allow cos 8 bi to be a free parameter. It is known from a comparison to 
the total Z width, however, that cos 6^ < 0.6 ^H]- We use the code Msmlib from G.Ganis 
to compute the sbottom cross sections [T7j, and fix the sbottom mass to 6 GeV. 

To examine the dependence of T on cos 6^ , we first fix a = 1. A shown below, this 
is close its optimal value. The results of the minimization of T are as one would expect. 
The first region is dominated by photon exchange, and so there is no sensitivity to cos 6 bi . 
The second retains a very mild sensitivity to cos9~ b , favoring cos 6~ b ~ 1, since the excess 
in this region is large. Region 3 is dominated by Z exchange, and so the dependence of the 
cross section on cos 6~ b is strong. It clearly favors a reduced coupling to the Z, and the limit 
cos6 b < 0.6 is easily obtained. Finally, the data from region 4 have a mild sensitivity to 
cos9 b , similar to what is seen for region 3. 

We sum the negative log-likelihood functions for the four regions and plot the sum as 
a function of cos 9 b , as shown in Fig. El This variation is dominated by the Z pole data. 
The preferred value is cos 6^ = 0.18, which gives a small but non-zero coupling to the Z, as 
required by the small excess in cx° ad . A zero coupling corresponds to cos% = 0.38, and is 
not excluded by the data - it is disfavored by only 1.3 standard deviations. 

We now fix cos9 b = 0.18, and return to the procedure to evaluate the significance of 
the excess as described in the last section. Specifically, we allow for an overall theoretical 
uncertainty, and we take correlations among measurements into account. In the present 
context we are comparing the quality of the description of the data for SM + a light sbottom 
to the SM alone. We vary a freely for the sum of negative log-likelihoods, and obtain 
dbest = 1.17 for which jF(a best ) = 50.83, to be compared to ^(0) = 60.03. In the latter case, 
the systematic uncertainty on the SM prediction allows an enhancement of the prediction 
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by 0.87%- i.e., Pbest = 1-0087 in Eq. ()12j) . This difference in T corresponds to 4.3 standard 
deviations, for a probability of 9.0 x 10~ 6 . 

It is important to note the following three points: 

1. All four regions contribute positively to AJF at c^best- 

2. The best overall point, «best, is not disfavored by any of them, and is statistically 
consistent with the best point for each individual region. 

3. The best overall point is consistent with a = 1. 

To illustrate these points graphically, we temporarily revert to the definitions of T which 
neglect theoretical uncertainties and correlations, and plot the variation of T for each of the 
four regions and their sum - see Fig. 

The minima of the four parabolas fall closer to a = 1 than to a = 0. This indicates 
there is a consistent excess in all four regions. The sum of the four parabolas has a sharper 
minimum at a = 1.2. The rise in the negative log-likelihood passing from a = 1.2 to a = 
corresponds, when interpreted as a number of standard deviations, to 6a. (Recall that, 
after taking theoretical uncertainties and correlations among measurements into account, 
this number is reduced to 4.3<r.) In this sense one can say that the data prefer a description 
which includes a light sbottom, with a high level of statistical significance. If the parabola 
were centered at a = instead, one would say that this model for a light sbottom is excluded. 
However, one cannot say that the data exclude the Standard Model. The contribution to 
e + e~ — > hadrons of SM processes is indisputable, while the contribution from light sbottoms 
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is speculative. The results of this analysis do not prove the existence of light sbottoms, rather, 
they suggest there may be new processes contributing to the hadronic final state. 

For illustration, we plot the difference between the rebinned measured values and the 
SM prediction as a function of y/s, as shown in Fig. |HJ The shaded band centered on zero 
indicates the conservative error we have assigned to the SM prediction. The heavy smooth 
curve indicates the cross section for a light sbottom (M^ = 6 GeV) when cos 6^ = 0.18. 
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Figure 7: Negative log-likelihood as functions of a. The mixing angle has been fixed by 
cos &i =0.18. The four lower parabolas show the contributions from each of the four regions, 
as marked. Their sum results in the uppermost parabola, which has a minimum close to 
a = 1, corresponding to the sum of SM and hi cross sections. This point is preferred over 
the SM alone a = 0. 
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Figure 8: Comparison of the data to the prediction of a light sbottom (Mi = 6 GeV) with 
cos8i = 0.18. The points represent the data, rebinned in the range 20 GeV < ^/s < 75 GeV. 
The upper plot shows the full range of ^Js, while the lower plot shows the Lep 2 data alone. 
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region 


«best 


-^(ttbest) 


^(0) 


S.D. 


1 


0.50 


21.89 


25.23 


2.58 


2 


0.27 


17.31 


37.50 


6.35 


3 


0.67 





1.404 


1.68 


4 


0.0007 


2.46 


6.12 


2.70 



Table 6: Results of the minimization of the negative log-likelihood function, T , as a function 
of a, for the case of a heavy Z' . 



5.2 A Heavy Z' Boson 

We consider the exchange of a third boson, call it Z', in the process e + e~ — > hadrons. 
For simplicity, we assume that there is no interference between this Z' and the photon and 
SM Z boson. If the mass of this Z', M z >, were much higher than the -^/s probed by these 
measurements, then cross section would rise approximately linearly in s: 

A = lna* QED C s («) 2 + (a'J 2 ) {{v' f f + {a' f f) . 

(13) 

The factor A contains the coupling and other factors. 

We use the expression Eq. © with y NP (s) = (10 pb) • (s/ (300 GeV) 2 ). A summary of the 
results for the four regions is given in Table El The combined likelihood gives a ne t = 0.00072, 
for a AjF = 3.74, corresponding to 2.8 standard deviations. This model makes only a small 
improvement over the SM. The rising s dependence pits the Lep 2 data against the data in 
regions 1 & 2 - one has either ?/ NP too high in region 4 or too small in regions 1 & 2. This 
model is distinctly less successful than the other two. 

6 Summary 

We have examined the measurements of a(e + e~ — *> hadrons) across a wide range of energies 
and experiments. None of these measurements shows a significant excess over the predictions 
of the Standard Model. Taken together, however, there is a clear trend, which has been 
quantified in a number of ways. 

Motivated by this apparent excess, we compared the data to a model for the production of 
light scalar 6-quarks, of the type suggested by the Argonne group [T5] . As a specific example, 
we have considered sbottoms with a mass of 6 GeV and a mixing of left- and right-states 
parametrized by cos6^ = 0.18. This mixing gives a relatively small but non-zero coupling 
to the Z boson. 

When all measurements are taken to be independent, and when theoretical uncertainties 
in the SM prediction are ignored, the light sbottom model is 'preferred' over the SM by 6.1cr. 
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s-Mf, + iM z > Y z , 
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When correlations among measurements are taken into account, this significance drops to 
5.2a. Finally, when a very conservative theoretical uncertainty is folded into the calculation, 
the resulting significance is 4.3(7, which corresponds to a probability of 9 x 10~ 6 . 

Alternative models are less successful in describing the data. A new neutral gauge boson 
is disfavored, and a cross section with a generic 1/s dependence does not explain the slight 
excess at Lep 1 (though it still provides an improvement over the SM with a significance of 
3.9(j). 

We conclude that the description of the data on e + e~ — > hadrons is significantly improved 
if Standard Model processes are augmented by the pair production of light scalar particles 
of charge ±1/3 which appear as hadronic final states. 
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